Observable estimation of entanglement of formation and quantum discord for 
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We present observable lower and upper bounds for the entanglement of formation (EOF) and quantum dis- 
cord (QD), which facilitates estimates of EOF and QD for arbitrary experimental unknown states in finite- 
dimensional bipartite systems. These bounds can be easily obtained by a few experimental measurements on a 
twofold copy g® g of the mixed states. Based on our results, we use the experimental measurement data of the 
real experiment given by Schmid et al. [Phys. Rev. Lett. 101, 260505 (2008)] to obtain the lower and upper 
bounds of EOF and QD for the experimental unknown state. 
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1. INTRODUCTION 

In recent years it has become more and more clear that 
quantum correlations, including quantum entanglement and 
quantum discord, are basic resources in quantum-information 
processing [1]. Therefore, the quantification of quantum cor- 
relations becomes fundamental problems in quantum infor- 
mation theory. However, quantum correlations are not yet 
fully understood. One of the most important entanglement 
measures is the entanglement of formation (EOF) [2]. For 
pure state \(f), it is defined by Ef{\(f)) - Sv{qa), where 
S vio) - -Tr(p log2 q) is the von Neumann entropy and qa = 
TiBi\(p}{f\) is the reduced density matrix of subsystem A. For 
mixed state g, the EOF is defined by the convex roof. 



Epie) = min V piEF{\ipi)), 



(1) 



where the minimum is taken over all possible ensemble real- 
izations Q = 2; P/IVfX^il with Pi > and 2; P/ = 1- 

The quantum discord (QD) provides a measure for the 
quantum correlation beyond entanglement [3,4]. It is believed 
that in certain quantum computing tasks there are still quan- 
tum advantages in the absence of entanglement. One typical 
example is the deterministic quantum computation with one 
qubit [5] in which the QD is proposed to be the reason for the 
quantum speedup [6]. Therefore, the QD has attracted much 
interest in quantum information theory [7-25]. The QD is de- 
fined to be [3,4], 

DAie) = min V pkS v(gB\k) + S v(gA) - S vig), (2) 

[Et) ~^ 
k 

where we perform the positive operator-valued measures 
(POVMs) [Ek] on subsystem A and pk - Tr(£'t i8> 1^), QB\k — 
T^AiEk ® lg)/pk- The minimum in Eq. (2) can also be taken 
over all the von Neumann measurements [3], and these two 
definitions coincide in the case of zero quantum discord. 



Interestingly, there is a close relationship between the EOF 
and QD, namely the Koashi- Winter relation [26]: 



l^AigAB) -Sv(gA) +Sv{gAB) = Ef(gBc), 



(3) 



* cqtzcj@nus.edu.sg 
^ phyohch@nus.edu.sg 



where Ef(gBc) is the EOF of ^bc and gsc = Tr^diAXfAUBc) 
with \4t)ABC being a purification of gAB- 

Although both the EOF and QD were introduced many 
years ago, they are notoriously difficult to compute because 
of the minimization. Until now, only a few analytical results 
for the EOF and QD have been carried out, such as EOF of 
two-qubit states [27, 28], isotropic states [29], Werner states 
[30], and QD of two-qubit Bell diagonal states [31], rank-2 
two-qubit states [32] and gaussian states [33]. In order to es- 
timate entanglement and QD, several lower and upper bounds 
of entanglement measures and QD have been proposed [34— 
43]. However, there are only a few bounds for EOF and still 
no general results of lower and upper bounds of QD for ar- 
bitrary finite-dimensional bipartite states. Therefore, general 
results of analytical lower and upper bounds for the EOF and 
QD are imminently needed. 

In this paper, we present observable lower and upper 
bounds for EOF and QD, which can allow estimates of EOF 
and QD for arbitrary experimental unknown states in finite- 
dimensional bipartite systems. These bounds can be easily 
obtained by a few experimental measurements on a twofold 
copy g®g of the mixed states. Based on our results we use the 
experimental measurement data of the real experiment given 
in Ref. [37] to obtain the lower and upper bounds of EOF and 
QD for the experiment unknown state. 



II. OBSERVABLE LOWER AND UPPER BOUNDS FOR 
ENTANGLEMENT OF FORMATION 



Before embarking on our main results, let us briefly review 
the von Neumann entropy, the linear entropy, and the lower 
and upper bounds of squared concurrence. The von Neumann 
entropy, generalized from the classical Shannon entropy to 
quantum states, is defined by 

S v{g) = -Tr(^ log2 e) = - ^ y"/ logi y"; = Hiji), (4) 



where yu, are eigenvalues of g and ft is the Schmidt vector 
{yUi,jU2, ■ ■ ■ ,//rf). The Hnear entropy, based on the purity of 
a quantum state, can be defined as 



Sdg) = 2(1 



V) 



2(1 -X^/)' 



(5) 



where for simplicity we have added a coefficient "2". The 
lower and upper bounds of squared concurrence for arbitrary 
finite-dimensional bipartite states introduced in Refs. [38, 39], 



Tr{g®gV,)<[C{g)]-<Ti(g®gKi) 



d(1) 



3(2) 



where V, = 4{P'y - P^') » P':', V: 

3(1) ^11(2) „„A I'. _ /i^-n(l) ^ d(2)\ ,„4fU d(') ,'d(') 



4P"> ® (P^i' 



(6) 



P?), 



Ki = 4P'i ' ® 1^2' and ^2 = 4(10 ® P^^O with P*!' (P';0 being 
the projector on the antisymmetric (symmetric) subspace of 
the two copies of the /th subsystem. It is worth noting that the 
lower and upper bounds can be expressed as follows. 



TT(g®gVi)-- 


-- 2(Tr£,2 - TTgl), 


(7a) 


TT(g®gV2)-- 


-- 2(TTg^ - TTgl), 


(7b) 


Trig®gKi)-- 


= 2(1-Trei), 


(7c) 


Tiig^gKi)-- 


--2(l-TTgl). 


(7d) 



Based on these bounds of the squared concurrence, we shall 
provide bounds for the EOF. 

For simplicity, we use the two denotations co{g) and ca{g). 
Here co{g) denotes the convex hull of the function g, which is 
the largest convex function that is bounded above by the given 
function g. Conversely, ca{g) denotes the smallest concave 
function that is bounded below by the given function g. The 
two denotations have been used to obtain explicit expressions 
and bounds for the EOF [29, 40^2]. 

Theorem 1. — For any m ® n {m > n) quantum state g, its 
entanglement of formation Ep(jg) satisfies 



co[Pf (A5)] < Epig) < ca[F[;'(A6)], 



(8) 



where A^ = max{0, 2(Tr£)2 - Trg^), 2(Tr£)2 - Txgl)), K(, = 
min{2(l - Tr£)^),2(l - Tr^)^)) (specially Ae = min{2(l - 
Trp^), 2(1- Trgi^), 1 - Tr^)^ - Tr^;^ + Txg^] for two-qubit states 
[39]) and 



Rf{X) = H2[ka{X)] + ka{X)\og2k, 



a{X) = [1 + ^]l-{k+\)A^I{2k)]l{k+ 1), 
^f(T) = H2[y{T)] + [1 - y(T)] log2(flf - 1), 

j{t) = [1 + yl{d-lY-d{d-\)Tl2]ld, 



(9) 
(10) 
(11) 
(12) 



with H2(x) - -X log2 x-{\- x) log2(l - x) being the standard 
binary entropy function, k - \_2I{2-A~)\ where [-J denotes the 
floor function, A e [0, V2(<i- l)/d] and t e [0, 2(rf - 1)/^]. 
Proof.- We first find the minimal admissible H{jl) for a 

given A - J2(l - 2,/U^), and the maximal admissible H(iT) 



for a given t = 2(1 - 2, fj.j). Consider the following two func- 
tions. 



R^^\A) = min Ih(P)\A = . 2(1 - ^/z^) 
{H{fT}\T = 2{l-Y,f^l)}, 



Fr,\r) 



maxi 

P 



(13) 



(14) 



As shown in Appendix, the minimal H(jj) versus A consists 
of d - I segments and the ^th segment corresponds to ft in 
the form {?,■■■ ,f,l - kt,0,--- ,0) for t e [l/(k + l),l//t], 
and the maximal H(jl) versus t corresponds to ft in the form 
{t,(l-t)/(d-l),--- , (1-0/(^-1)1 for fe [l/t/,1]. Therefore, 
one can obtain Eqs. (9) and (11). 

Suppose that we have already found an optimal decompo- 
sition Yjj Pj\^j)(^j\ for 6 to achieve the infimum of Ef(g), 
then Efig) = Tij PjEpilt/'j)) by definition. Since co[/?^ (A)] 
is a monotonously increasing convex function and satisfies 
co[P^ (A)] < H(ft) for a given A, one thus has 



:co[4"'(2/V')]>co[P<;"(A5)], 



(15) 



where we have used 2y/?y/l-' > A5 which has been proved in 
Ref. [38]. Meanwhile, since ca[F~^ (r)] is a monotonously 
increasing concave function and satisfies ca[F^ (t)] > H(ft) 
for a given t, one thus has 

Epig) = 2/^,£f(|^,)) = 2/7,//(^^) < 2/7,ca[F<;'(W)] 



ca[F*;>(2/7,r>)]<ca[P;;>(A6)], 



(16) 



where we have used Y^jPF' ^ ^6 which has been proved in 
Ref. [39]. ■ 

Remark 1. Actually, Ref. [42] has got a similar result. 



which is equivalent to co[P'^ (A^)] < Ep{g) < ca[Fy (Ae)] 



where F^ (r) is defined in the appendix. However, our 
lower bound co[P^''(A5)] is better than their lower bound 
co[F£' (Aj)]. Furthermore, we will give analytical results of 
co[P^ (A5)] and ca[pj (A5)] in the following, and our ana- 
lytical results are different from the results shown in [42]. 

Remark 2. The above theorem provides explicit lower and 
upper bounds of the EOF. Since A5, Ae (as well as Ai, ■ ■ ■ , A4 
shown in the next section) can be directly measured in exper- 
iments using two copies of the state, one can easily obtain the 
lower and upper bounds of the EOF without quantum state 
tomography. 

As introduced above, co{g) is the largest convex function 
that is bounded above by the given function g, and ca(g) de- 
notes the smallest concave function that is bounded below by 
the given function g. From Eq. (9) one can obtain the explicit 
expression ofco[Rf^ '(/I)] which also consists ofd-\ segments, 
meanwhile from Eq. (11) the explicit expression of ca[F[^ (t)] 

can be obtained. Therefore, co[7?^ (A)] and ca[F[^ (t)] read as 
follows (the detailed calculations have been given in the ap- 
pendix). 



co[<'(^)] 



ca[F[f (t)] = 



f H^{^^^), 



Ae [0,1], 



log2(t+l)-log 



m[^- ^y^o^.ik), Ae[l,^l 



H2[y{T)] + [l-y(T)]\og^(d-\), 



T€ [0, 



Ad-b 1 
rf(i/-l)J' 



^l0g2(t/-l)- 



d(d-2) 



^ log2((i - 1) + log, (i, T € [ 



4t/-6 2(t;-l) -| 
d(d-l)' d J 



(17) 



(18) 



where j{t) 
A^)\. 



[1+ V(^r^^3^^^^^(crn>72]/Jand/t = L2/(2- 



III. OBSERVABLE LOWER AND UPPER BOUNDS FOR 
QUANTUM DISCORD 

Actually, from the proof of Theorem 1, we can also obtain 
the lower and upper bounds for QD. 

Theorem 2. — For any m®n quantum state g, the quantum 
discord satisfies 

OAig) > RfiAi) - r';;"'\A2) + co[4"'(A3)], (19) 



where A? 



2(1 - Txel\ A? 



.1 - ^v. ^ ..^^y, -2 - 2(1 - Txe\ A^ = 
max{0, 2{Tiq\ - Trp|), 2(Tvq\ - Txg^)] and 

Rf{A) = H2\Ji{A)] + [1 - piA)] log2(flf - 1), (20) 

p{A) = [1 + yl{d-l)^-d{,d-l)A^I2\ld, (21) 

with /I e [0, ^J2{d- I) Id]. 

Proof. — We first find the maximal admissible H{fj) for a 

given A - J2(l - 2,;U^). Consider the following function. 



R^^Ha) = max \ H{f})\A = 



(I 



(22) 



As shown in Appendix, the maximal H{jT) versus A corre- 
sponds to ft in the form (f, (1 - t)/{d - 1), ■ ■ ■ , (1 - t)/(d - 1)) 
for t e [l/d,l]. Therefore, one can obtain Eq. (20). 

In order to find the lower bound of the minimization term 
in DAigAs), we use the Koashi-Winter relation Eq. (3). Sup- 
pose that we have already found an optimal decomposition 
YjjPi\<pj){if>j\ for gsc to achieve the infimum of Ef(gBc), then 
Ef{Qbc) = I,jPjEf{\iI>j}) by definition. Since co[R'^1^\A)] 
is a monotonously increasing convex function and satisfies 
co[Rj^ (A)] < H(Jl) for a given A, one thus has 

i J ;■ 

> co[7;*;"(^ pjA^)] > co[7?["'(A3)], (23) 

;■ 

where we have used YijPj^^ ^ A3 which has been proved in 
Ref. [38]. Moreover, from the definitions in Eqs. (13) and 
(22), we can obtain that 



Sv(gA)>Rf\Ax\ 
SvigAB)<R^r^(A2l 



(24) 
(25) 



I 

which together with Eq. (23) gives exactly Eq. (19). ■ 

In the following, we shall also present an explicit upper 

bound for the QD, which is observable as well as the lower 

bound shown in Theorem 2. 

Theorem 3. — For any m ® n quantum state g, the quantum 

discord satisfies 

i)A{Q) < R^u\Ai) - Rf'"\A2) + ca[F[;\A4)], (26) 

where A4 = min{2(l - Trg^), 2(1 - Tvg-)}. 

Proof. — Similar to Theorem 2, we also use the Koashi- 
Winter relation Eq. (3). Suppose that we have already found 
an optimal decomposition Y^j Pj\4'j){4'j\ for gsc to achieve the 
infimum of EfigBc), then Efiggc) = Yjj PjEf{\^j)) by defini- 
tion. Since ca[Fy (r)] is a monotonously increasing concave 
function and satisfies ca[Fy (r)] > //(/?) for a given t, one 
thus has 

Ef(Qbc) = Xi PjEf(\4>j)) = Y, Pj^^i^'^ - Z PiCa[4"^(T0] 

j J j 

< ca[F^;;\J] PjT^)] < ca[4"'(A4)], (27) 

j 

where we have used YijPj'^^ ^ A4 which has been proved in 
Ref. [39]. Moreover, from the definitions in Eqs. (13) and 
(22), we can obtain that 



Sv{QA)<Rf{A,), 

^(mn). 



(28) 
(29) 



Sv(gAB)>Rr(^2), 

which together with Eq. (27) gives exactly Eq. (26). ■ 

Remark 3. The lower and upper bounds of QD are also 
valid for the minimum in Eq. (2) being taken over all the 
von Neumann measurements, since the QD using POVMs is 
always smaller than or equal to the one using von Neumann 
measurements, and SjPjT-' < A4 used in Eq. (27) still holds 
for the case using von Neumann measurements. 



IV. EXAMPLES 

In this section, we will present several examples using our 
lower and upper bounds of the EOF and QD. 

Example 1. — Considering pure states in 2 ® n systems, one 
can observe that our lower and upper bounds of EOF and QD 
coincide. Because for these pure states we have 

/?f (Ai) = ;?g*(Ai) = co[Rf{A,)] = c^[Ff{Ae)l (30) 



R 



(2n)/ 



yiMt 



vj, XA2) = /?r (^2) = 0, 

co[/?f (A3)] = ca[F*;"(A4)] = 0, 



(31) 
(32) 




0,0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 

Observable lower bound of EOF 

FIG. 1. Observable upper bound of the EOF versus its observable 
lower bound for (2 ® 4)-dimensional random states with different de- 
grees of mixing: (a) shows weakly mixed states (0.10 < yi-Trp^ < 
0.11), (b) displays intermediate mixing (0.20 < -^1 - Trp^ < 0.21), 
and (c) corresponds to strongly mixed states (0.30 < -y'l - Tig- < 
0.31). The dashed lines denote the lower bound. 



it can be directly obtained that 

EAg) = Da{0) = <*(Ai) = Rf(Ai). (33) 

Therefore, for the pure sates in 2 18>« systems the EOF and QD 
can be directly measured using our lower and upper bounds. 

Example 2. — Let us simulate the lower and upper bounds 
of the EOF and QD for mixed random states of (2 i8> 4)- 
dimensional systems. The mixed random states are generated 
with different degrees of mixing. Using our theorems, the up- 
per bounds versus the lower bounds of EOF and QD are de- 
picted in Fig. 1 and Fig. 2, respectively. For weakly mixed 
states, the lower and upper bounds provide an excellent esti- 
mation of EOF and QD; for strongly mixed states, the bounds 
also provide a region for EOF and QD. 

Example 3. — Using the experimental measurement data of 
the real experiment shown in Ref. [37], one can obtain the 
lower and upper bounds of the EOF and QD for this experi- 
mental unknown two-qubit state g^^p. Ref. [37] has obtained 
the following measurement results. 



P- 



0.208, P-+ = 0.050, SJ+ 



0.061, 



(34) 



where p = {P[_' i8> P ), p-+ and p+^ being defined simi- 
larly. Using Eqs. (7) and (34), one can calculate that 



Ai 

A3 

As 



1.016, 
0.210, 
0.795, 



A2 
A4 

Afi 



0.666, 
0.444, 
0.832. 



(35) 



Therefore, the lower and upper bounds of the EOF and QD 
for this experimental unknown state are 



0.715 <£F(e«p)< 0.875, 
QAlQ<DA{Qexp)< 1-023. 



(36) 
(37) 




0.5 1.0 0.0 0.5 1.0 0.0 

Observable lower bound of QD 



1.0 



FIG. 2. Observable upper bound of the QD versus its observable 
lower bound for the same states used in Fig. 1. 



Although the upper bound in Eq. (37) gives an invalid bound 
[44] , the lower bound provides an efficient lower bound of QD 
for the experiment unknown state. The reason for the invalid 
upper bound comes from the experimental error. It is worth 
noticing that Ai = 2(1 - Trp^) which should be less than or 
equal to 1 for two-qubit states. However, in this experiment 
one has Ai = 1.016 because of the experimental error. There- 
fore, the real upper bound of QD without the experimental 
error is probably close to and smaller than 1. 



V. DISCUSSION AND CONCLUSION 

Let us discuss the realization of the measurements for our 
lower and upper bounds in real experiments. As shown in 
Ref. [37], the antisymmetric projector P for two-qubit states 
takes the particularly simple form P - \il/'){\l/^\ with \^~) — 
(|01) - |10))/ "V^ being the singlet state. It can be measured 
by several methods using two copies of g in photonic sys- 
tem. For example, Huang et al. [45] have used a Hong-Ou- 
Mandel interferometer [46] to project two photons onto the 
singlet state, and Refs. [37, 47] have used a controlIed-NOT 
gate to distinguish the Bell states, since the controUed-NOT 
gate can transform the Bell states to separable states and dis- 
tinguishing separable states is much easier than distinguishing 
the Bell states. 

In this work, we have presented observable lower and up- 
per bounds of the EOF and QD. These bounds can be used 
to estimate EOF and QD for arbitrary experimental unknown 
states in finite-dimensional bipartite systems. One can easily 
obtain these bounds by a few experimental measurements on 
a twofold copy g® g of the mixed states. Furthermore, based 
on our results we have used the experimental measurement 
data of the real experiment given in Ref. [37] to obtain the 
lower and upper bounds of EOF and QD for the experiment 
unknown state. 
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APPENDIX 



In this appendix, we exhibit the details to obtain the expres- 
sions of F[f*(T), Rf{A), Rf(A), co[R'^^\A)] and ca[F[f'(T)]. 
The main idea is to illustrate the difference between von Neu- 
mann entropy and the linear entropy (or the square root of the 
linear entropy). 



A. Calculation of f'^ (t), R'f^ (A), i?'^' (A) 



In the following, we shall seek the highest and lowest von 
Neumann entropies consistent with a given value of linear en- 
tropy. Actually, this problem is equivalent to seek the max- 
imal and minimal H(J}) = - Zi/A'/ log2jU, for a given r = 
2(1 - Yiil^J)^ ^nd the later one is a classical problem which has 
been solved in Refs. [48, 49]. We use F^{^\t) and ^^'(t) to 
denote the maximal and minimal H{fj) for a given r (although 
Fj (t) has not been used in the main text, we introduce it for 



the sake of completeness), i.e.. 



Ff{T) = max {h{jD\t = 2(1 - V ^iH, (38) 

d 
Ff{T) = mm [Hift)\T - 2(1 - J]/zf)}, (39) 



(=1 



where d denotes the dimension of the Hilbert space. 

As shown in Refs. [48, 49], the maximal H(jt) versus t 
corresponds to fl in the form 



1 



^i 



d-l 



i-4) forf€[i 1] (40) 

d-l d 



with d- 1 copies of(l-t)/{d- 1) and one copy off. Therefore, 
the maximal //(/?) and corresponding t are 



//(f) = -flog2f-(l-f)log2 



1 



d-l' 



T(t) = 2[1 - f2 - 



d-l 



(41) 
(42) 



In order to show the maximal H(p) versus t, we need the in- 
verse function of T{t). After some algebra, one can obtain that 



Kt) = - + 
d 



1 (d-\f (d-\)T 



d'~ 



2d 



(43) 



with re [Q,2{d- I) Id]. Substituting Eq. (43) into Eq. (41), 
we can get the expression for Fj^(t), i.e., Eq. ( 1 1 ) in the main 
text. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.1 



FIG. 3. (Color online) The linear entropy versus the von Neumann 
entropy. 50000 dots represent randomly generated states with d = A. 
Rank-2, rank-3, rank-4 states are denoted by the blue, red, black dots, 
respectively. The upper boundary is a smooth curve, whereas the 
lower boundary consists of three segments. 



The minimal H{}t) versus t consists of li - 1 segments and 
the feth segment corresponds to ft in the form [48, 49] 



fl-{t,--- ,t,l -kt,Q, 



,0) forfe[-l-,|] (44) 
k + 1 k 



with k copies of t, d-k- I copies of and one copy of 1 - kt. 
Therefore, the minimal H(ft) and corresponding r are 

H{t) = -kt \og2 t-il-kt) log2(l - kt), (45) 

T(t)^2[l-kt^-{l-ktf]. (46) 

In order to show the minimal H{jj) versus t, we also need the 
inverse function of T{t). After some algebra, one can obtain 
that 



t(T) 



1 



1 



k+1 k+l 



{k+l)T 



2k 



(47) 



with r e [0,2(d - l)/d]. Substituting Eq. (47) into Eq. (45), 
we can get the expression for F^ (t), i.e., 

^f(r) = //2[M(T)]+M(T)log2^, 



6iT) = [1 + yll-ik+l)T/(2k)]/{k+ 1), 



(48) 
(49) 



with/t ^ !,■■■ ,d- I. Since f e [l/(k+ l),l/k] for the Wh 
segment, we have t € [2{k - I) Ik, 2k/(k -Hi)]. Therefore, one 
can obtain k - [2/(2 -t)J where [-J denotes the floor function, 
i.e., [x] is the largest integer not greater than x. 

We can simulate the lower and upper boundary of the re- 
gion in Sl versus Sv plane. 50000 dots for randomly gen- 
erated states with d = 4 are displayed in Fig. 3 (a similar 
figure has also been shown in Ref. [49]). The upper boundary 
corresponds to the function F^ (t), and the lower boundary 
corresponds to F^ (t) with three segments. 

Similarly, we use R]j (A) and Ri^ (A) to denote the max- 
imal and minimal H(jT) - -Yjif^i^'^Ei/^i for ^ given A - 



w 




B. Calculation of co[iff (^)] 



FIG. 4. (Color online) The square root of linear entropy versus tlie 
von Neumann entropy based on the same states used in Fig. 3. 



y(d). 



We have already obtained the expression for R^ (A) in the 
above subsection, which comprises d - I segments with k = 
I,--- ,d-\. In order to obtain co[7?^ (/I)], we need to find the 
largest convex function that is bounded above by R^ (A). It is 
worth noting that the first segment of /?^ (A) (i.e. ^ = 1 in Eq. 
(9)) is convex, whereas the rest segments (k - 2, ■ ■ ■ ,d- 1) 
are concave functions, and one can take Fig. 4 as an exam- 
ple. Therefore, the first segment of co[/?^ (A)] is R^ (A) with 

k = 1 itself, i.e., H2(l/2 + Vl - A^/2) for A e [0, 1]. The kth 
segment of co[R['(A)] (k - 2,- ■ ■ ,d - I) is the line between 
two points: ( V2(fc- l)/k, logj k) and ( yJ2k/(k+ 1), log2(A: + 
1 )), i.e., {A - V2(^-l)/^)log2(l + lA)/(V 2fc/(^+l) - 
V2(/t- V)lk) + logj/t for /I e [1, ^2(d- V)ld\. Therefore, 
the explicit expression of co[7?^ (/I)] is Eq. (17). 



^ 



2(1 - Yui^^i)^ i-S-, Eq. (22) and Eq. (13) in the main text. 
Since A is the square root of t, the forms of /? for the maxi- 
mal and minimal //(/T) versus A are the same as the forms for 
the maximal and minimal //(/7) versus r, i.e., Eq. (40) and 
Eq. (44), respectively. Therefore, the maximal //(/I) and cor- 
responding A are 



i/(f) = -flog2f-(l-01og2 



1 -t 
d-V 



Aii)^ ,/2[l-f2-^L4!]. 



d-\ 



and the inverse function of A{t) is 



(50) 
(51) 



C. Calculation of ca[F'f'(T)] 



M)/ 



We have already obtained the expression for F^ (t) in 
the above, which is a smooth function. In order to obtain 
ca[Fy'(T)], we need to find the smallest concave function that 
is bounded below by F^ (t), which is similar to the situation 
inRef. [41]. 

When d -2, since F^ (t) is a concave function, ca[F^ (t)] 
is F^ (t) itself. When d>3,we first prove that there is one 
and only one point To between and 2(d - l)/d such that 



7W' 



7('0/ 



' u (''"o) - 0- The second derivative of Fjj (t) with respect 
to r is 



^ ' d y d^ 2d 



(52) 



with A € [0, V2((i- l)/d]. Substituting Eq. (52) into Eq. 
(50), we can get the expression for R^ (A), i.e., Eq. (20) in the 
main text. Furthermore, the minimal H(jT) and corresponding 
A are 



H{t) = -kt log2 1 -{1- kt) log2( 1 - kt), 
A(t) = ^J2[l-kfi-(l -kt)^, 

and the inverse function of A{t) is 



f(^) = 



k+l k+l 



(k + 1)A^ 



2k 



(53) 
(54) 



(55) 



with A 6 [0, yJ2(d- l)/d]. Substituting Eq. (55) into Eq. 
(53), we can get the expression for R^ (t), i.e., Eq. (9) in the 
main text. 

Fig. 4 shows the simulation for the lower and upper bound- 
ary of the region in V^T versus S y plane. The same randomly 
generated states has been used as Fig. 3. The upper boundary 
corresponds to the function R^ (A), and the lower boundary 
corresponds to R, (A) with three segments. 



7(dy 



(t) 



■log 



(d-l)7(T) 
_ dyir) - 1 1 - -/(t) 

x(y'(T))' log2 e. 



1 



r(T)(i-r(T)) 



(56) 



where Y{t) - -{d - l)/(4<iy(r) - 4) and Y'{t) = 



^d(r'(T))'/{d ■ 



(0) = 



lim. 



7W' 



1). From Eq. (56) we can obtain F\j 

oF[f(2(d- 



'^.Ff/\2{d-\)ld) 



lim. 



l)/d - e) - +00 and tq is the solution of /(r) = g(T) where 



fir) = log 



id-l)y{T) 
l-r(T) ' 
dyJT) - 1 
dr(T)(l-y(T)y 



(57) 
(58) 



For simplicity we define y = 7(r)/(l-y(T)) withy € (l/(d- 
1), +00) because of t e (0, 2(d - l)/d). Thus, one can directly 
consider the functions / and g with respect to y. 



f(y)^log[{d-l)yl 

1 2 1 

8iy)-(i-::)y + ii-::)-:r- 

d d dy 



(59) 
(60) 



Therefore, tq corresponds to the solution of A{y) = f{y) - 




FIG. 5. (Color online) A typical figure of A(y) versus y when d = 3. 
In this case y, = 1/2, y2 = 1, ^3 = 2/3, and yo ~ 1.40845. 



g(y) - 0. It is worth noting that 

AO;) = log[{d - l)y] - (1 - i)y - (1 - ^) + 1, (61) 

a a ay 



1 1 1 
dy^ y a 



A"(y) = 



2-dy 
dy^ 



(62) 
(63) 



One can directly obtain that A'{y) = has two roots yi 



l/(d-l) and ^2 = 1 , and A" (y) = has only one root ^3 -2/d. 
Moreover, we have A(yi) = 0, A{y2) > for d > 3, and 
A'iy) > when y G (l/((i - 1), 1). Therefore, A(3') is a convex 
function when y e (l/((i - 1), 2/d), but A(y) is concave when 
y e [2/ii, +00). Thus, when d > 3 there is only one solution yo 
(corresponding to to) for A(y) = when y e {l/{d - 1), +00). 
A typical figure of A(y) has been shown in Fig. 5 when d - 3. 
The above analysis indicates that for d > 3 F^ (t) is con- 
cave (i.e. F^j (r) < 0) near t = but convex (i.e. F^j (t) > 

0) near r = 2(d — l)/d. In order to obtain ca[Fy (t)], we 
should solve the following equations: Let 



£'line(T) = a 



2(d- 1) 



log2 d 



(64) 



be the line crossing through the point (2{d - l)/d, log2 d). We 
solve 



7(d)( 



£line(T) = F'^'iT) 
d£'line(T") 



dF^^\T) 



dT 



dT 



— a. 



(65) 
(66) 



for a and r. The solution to the equations is unique: t - {4d - 
6)/{d{d-l))mda = ((i-l)/(2£/-4)log2(ii-l). Therefore, for 
^ > 3 ca[F<f (t)] is Fifir) itself when r G [0, (4-d - 6) / (did - 



1))] and the line Eiine when r € [{4d-6)/{did-l)),2id- 
respectively. 
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